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We discuss possible origins of transverse spin asymmetries in hadron-hadron collisions and 
propose an explanation in terms of a chiral-odd T-odd distribution function with intrinsic 
transverse momentum dependence, which would signal a correlation between the transverse 
p?«.' [ spin and the transverse momentum of quarks inside an unpolarized hadron. We will argue 

that despite its conceptual problems, it can account for single spin asymmetries, for example in 
pp^ —>■ n X, and at the same time for the large cos2(/) asymmetry in the unpolarized Drell-Yan 
cross section, which still lacks understanding. We use the latter asymmetry to arrive at a crude 
model for this function and show explicitly how it relates unpolarized and polarized observables 
in the Drell-Yan process, as could be measured with the proton- proton collisions at RHIC. 
Cn , Moreover, it would provide an alternative method of accessing the transversity distribution 

^N ■ function hi. For future reference we also list the complete set of azimuthal asymmetries in 

the unpolarized and polarized Drell-Yan process at leading order involving T-odd distribution 
functions with intrinsic transverse momentum dependence. 









cd 



13.88.-fe; 13.85.Qk 



(N 
> 

in 

vn 

(N 

o 

Q>^ , I. INTRODUCTION 

r^ ■ Large single transverse spin asymmetries have been observed experimentally in the process pp^ —> tt X |l[| 

^ i' and many theoretical studies have been devoted to explain the possible origin(s) of such asymmetries. However, 
JL , one experiment cannot reveal the origin(s) conclusively and one needs comparison to other experiments. In this 
^t*' article we will use additional experimental results to propose an explanation in terms of a chiral-odd T-odd 
r^ , distribution function with intrinsic transverse momentum dependence and we contrast it to the more standard 
I!" ' theoretical proposals (cf. E|). In addition, it can account for the large cos2(/> asymmetry in the unpolarized 
.J^ , Drell-Yan cross section P,i|, which still lacks understanding. 

S^ ■ Unlike the chiral-even T-odd distribution function with intrinsic transverse momentum dependence as in- 

%^ I vestigated by [p|q], which depends on the polarization of the parent hadron, the chiral-odd function signals a 
correlation between the transverse spin and the transverse momentum of quarks inside an unpolarized hadron. 
But one can use the polarization of one hadron to become sensitive to the polarization of quarks in another, 
unpolarized hadron. In this way it could provide a new way of measuring the transversity distribution function 
hi. We propose two measurements that could be done at RHIC using polarized proton- proton collisions to 
study such a mechanism and to try to obtain information on hi. 

Apart from discussing the advantages of this proposal, we will discuss the theoretical difficulties connected 
to such a function. The function is actually the distribution function analog of the fragmentation function 
associated with the Collins effect [Q. Unlike the fragmentation function the distribution function is expected 
to be zero due to time reversal symmetry, unless one assumes some nonstandard mechanism to generate such 
a function, like for instance factorization breaking, which implies non-universality, or effects due to the finite 
size of a hadron, in which case one has the problem of systematically taking into account such effects in hard 
scattering factorization. 

The article is organized as follows. In Section II we discuss possible origins for transverse spin asymmetries. 
In Section III we elaborate on transverse momentum dependent distribution functions. In Section IV and also 
in the appendix, we give results for the leading order Drell-Yan cross section in terms of the T-odd distribution 
functions, for completeness taking into account contributions from Z-bosons. In Section V we discuss how one 
particular function can not only explain (in principle) the single spin asymmetries in the process pp^ — *■ tt X, but 
also explain the azimuthal cos 2(j) dependence of the unpolarized Drell-Yan cross section data ^j^ . In Section 
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VI we propose the measurements that could be performed at RHIC and that might uncover such an underlying 
mechanism. In Section VII we discuss the conceptual and theoretical problems related to T-odd distribution 
functions. 



II. ORIGINS OF TRANSVERSE SPIN ASYMMETRIES 

We will discuss possible origins of transverse spin asymmetries in the context of the following hard scattering 
processes, that either have been or will be performed. First we will go into the details of the single and double 
polarized Drell-Yan process Hi H2 -^ i£X, for which there is no data available yet. Then we focus on the 
single polarized process pp^ -^ tt X for which large single transverse spin asymmetries have been observed R]. 
We will also make use of knowledge of the unpolarized processes: n^ N -^ /i+ /i^ X and e+ e^ -^ hi ft,2 X. 



A. The polarized Drell-Yan process 



Transverse spin asymmetries in hadron-hadron collisions require an explanation that involves quarks and 
gluons. A large scale (the center of mass energy or the large lepton pair mass) allows for a factorization of 
such a process into parts describing the soft physics convoluted with an elementary cross section. The parts 
parametrizing the soft physics cannot be calculated within perturbative QCD. Let us first focus on the Drell-Yan 
process, i.e. lepton pair production in hadron-hadron collisions. 

In lowest order, i.e. the parton model approximation, the Drell-Yan process consists of two soft parts (in 
Fig. 1^ the leading order diagram is depicted g) and one of the soft parts is described by the quark correlation 
function $(Pi, Si;p) and the other soft part by the antiquark correlation function, denoted by ^{P2, S2; k): 
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FIG. 1. The leading order contribution to the Drell-Yan process 

As will be discussed in the next section one can decompose the quark momenta p and k into parts that are 
along the direction of the parent hadron, the so-called lightcone momentum fractions, and deviations from that 
direction. In case one integrates over the transverse momentum of the lepton pair one only has to consider the 
correlation functions as functions of the lightcone momentum fractions. 

The most general parametrization of the correlation function $ as a function of the lightcone momentum 
fraction x, that is in accordance with the required symmetries (hcrmiticity, parity, time reversal), is given by: 
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H^) = 7T ifii^)P + 9i{x) A75f + hi{x) Ts^t/'] 



(3) 



Other common notation is q for /i, Ag for gi and 5q or At<Z for hi. 

At this parton level one finds the well-known double transverse spin asymmetry 



Att oc \Sij,\ \S2t\ cos{(f>Si +(t>S2) hi{xi)hi{x2), 



(4) 



which has not yet been experimentally observed, but is one of the objectives of the polarized proton-proton 
scattering program to be performed at RHIC. This asymmetry is one of the possible ways to get information 
on the transversity distribution function hi. 

At the parton model level there are no single transverse spin asymmetries, but these might arise from cor- 
rections to this lowest order diagram. The corrections are of two types: the perturbative and the higher twist 
corrections. The first type depends logarithmically on the hard scale and the second type behaves as inverse 
powers of the hard scale. 

Assuming that single spin asymmetries arise due to perturbative contributions is conceptually the simplest 
option, since it assumes that the asymmetries actually occur at the quark-gluon level, i.e. they arise from ele- 
mentary subprocesses, and that going to the hadron level just involves convoluting the elementary asymmetry 
with (polarized) parton distributions. Typically this will yield single transverse spin asymmetries of the or- 
der asiTiq/y/s M which is expected to be sniallH. The perturbative corrections to the double transverse spin 
asymmetry Eq. (3) have been calculated in |0] and using the assumption that at low energies the transversity 
distribution function hi equals the helicity distribution function gi, it has been shown in Ref. |12[ that Att is 
expected to be of the order of a percent at RHIC energies. We will view this as an indication that perturbative 
QCD contributions are most likely not the (main) origin of large transverse spin asymmetries. 

Dynamical higher twist corrections to the parton model require expanding the correlation function $(a;) to 
include contributions proportional to the hadronic scale (typically the hadron mass), since these will show up 
in the cross section suppressed by 1/Q, where Q is a hard scale. At leading order in a^, i.e. (as)'^, but at the 
order 1/Q, one finds ||l^,0 no single or double transverse spin asymmetries^. 

Hence, in order to produce a large single transverse spin asymmetry one needs some conceptually nontrivial 
mechanism, since regular perturbative and higher twist contributions appear to be either small or absent. Two 
such nontrivial mechanisms are the soft gluon/fermion poles suggested by Qiu and Sterman mm and so-called 
time-reversal (T) odd distribution functions (cf. e.g. (l7J). Both of these mechanisms could produce a single 
transverse spin asymmetry At at order 1/Q. Recently it has been shown |lj] that their effects are identical in 
the Drell-Yan process, so in order to discriminate between them one must use other experiments as well. This 
asymmetry At has been estimated to be of the order of a percent at HERA energies (820 GeV, fixed target) 
|l8| . Let us remark that T-odd functions need not signal actual time reversal symmetry violation. We will 
discuss other options more extensively in the discussion at the end. 

In order to arrive at a single transverse spin asymmetry that is not suppressed by inverse powers of the hard 
scale, one can consider cross sections differential in the transverse momentum of the lepton pair. In that case 
one is sensitive to the transverse momentum of the quarks directly and in case this concerns intrinsic transverse 
momentum of the quarks inside a hadron, the effects need not be suppressed by 1/Q. The point is that if the 
transverse momentum of the lepton pair is produced by perturbative QCD corrections, each factor of transverse 
momentum has to be accompanied by the inverse scale in the elementary hard scattering subprocess, that is 
by 1/Q- But in case of intrinsic transverse momentum the relevant scale is not Q, but the hadronic scale, say 
the mass of the hadron. In processes with two (or more) soft parts, like the Drell-Yan process, the intrinsic 
transverse momentum of one soft part is linked to that of the other soft part resulting in effects, e.g. azimuthal 
asymmetries, not suppressed by 1/Q. These effects will show up at relatively low (including nonperturbative) 
values of Qt, where Qt = Or s-nd qt is the transverse momentum of the lepton pair. Studying the dependence 
of asymmetries on transverse momentum is another way to try to discriminate between possible origins for 
asymmetries. 

Returning to the parton model diagram and including intrinsic transverse momentum dependence in this pic- 
ture, one observes the following points. The effects will only show up if Qt is observed (i.e. not integrated over). 
If only T-even structures are included, several double spin azimuthal asymmetries are obtained, but no single 
spin asymmetries H]. So again one needs to include some nontrivial mechanism. Gluonic and fermionic poles 
have as yet not been considered with transverse momentum dependence (other than perturbatively produced), 
but would in any case appear in the cross section suppressed by a factor of 1/Q. However, the leading twist 
T-odd distribution functions with intrinsic transverse momentum dependence do yield single spin azimuthal 



^ Heavy quarks will appear at higher orders in as and have been shown to give rise to only small contributions (even) 
to the unpolarized Drell-Yan cross section jlO{ . 

^There is however a double spin asymmetry Alt, which involves one longitudinally and one transversely polarized 
hadron ||13l. A recent estimate of Alt using the bag model indicates that it is an order of magnitude smaller than the 
leading order asymmetry Att |15|]. 



asymmetries. We will be mainly focusing on the effects of such functions from now on. 



B. Pion production in pp^ scattering 

The large single transverse spin asymmetries that have been observed in the process pp^ ^ ttX |l|] require 
as said an explanation that involves quarks and gluons. Again one needs large scales (in this case also a large 
transverse momentum of the pion) to allow for a factorization of this process into parts describing the soft 
physics convoluted with an elementary cross section. For example, one contribution is coming from the diagram 
depicted in Fig. |[ 




FIG. 2. A contribution to the process pp^ -^ nX 

Assuming (as argued above) that perturbative and higher twist corrections (gluonic and fermionic pole con- 
tributions to this process have recently been investigated in Ref. |jl3]) are too small to generate the observed, 
large single transverse spin asymmetries, we will restrict ourselves to the transverse momentum dependent T- 
odd functions, in this case both distribution and fragmentation functions. The so-called Sivers H and Collins 
M effects are examples of transverse momentum dependent T-odd distribution and fragmentation functions, 
respectively. Like the transverse momentum of the lepton pair in the Drell-Yan process, the transverse momen- 
tum of the pion now originates from the intrinsic transverse momentum of the initial partons in addition to 
transverse momentum perturbatively generated by radiating off some additional parton(s) in the final state. 

Anselmino et al. m have investigated both the Sivers and the Collins effect as possible origins for the asym- 
metries as observed in Ref. [y. Both effects can be used to fit the data, which then can be tested using other 
observables. However, there are indications po[ from analyzing a particular angular dependence (a cos 2ip de- 
pendence |2l)) in the unpolarized process e'^ e^ ^ Z^ ^ mr X , where the pions belong to opposite jets, that 
the Collins effect, is in fact at most a few percent of the magnitude of the ordinary unpolarized fragmentation 
function. Therefore, it seems unlikely that the Collins effect is the main source of the single spin asymmetries 
of the pp^ -^ Ti X process. 

One other possible T-odd function that could be the source of the single spin asymmetries is the chiral-odd 
function hi , the distribution function analog of the Collins effect. It will be discussed extensively below for the 
case of the Drell-Yan process, but it can equally well be the source of single spin asymmetries in pp^ —^ tt X. 
On the other hand, the Collins effect itself will not contribute to the Drell-Yan process. 



C. The unpolarized Drell-Yan process 



The unpolarized cross section as measured for the process tt 
tungsten, using a 7r~-beam with energy of 140, 194, 286 GeV | 



N — > 11^ fi X , where N is either deuterium or 
I and 252 GeV [|, 
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shows remarkably large values of v. It has been shown [py22|] that its magnitude cannot be explained by leading 
and next-to-leading orderperturbative QCD corrections. A number of explanations have been put forward, like 
a higher twist effect [ESpj], which is the l/Q^ term discussed by Berger and Brodsky |25|. In Ref. |23] the 



higher twist effect is modeled using a pion distribution amphtude and it seems to fall short in explaining the 
large values as found for v. This higher twist effect would not be related to single spin asymmetries. 

In Ref. p2| factorization breaking correlations between the incoming quarks are assumed and modeled in 
order to account for the large cos 20 dependence. We will return to that extensively in section V. Another 
approach is put forward in Ref. pq] using coherent states. This can describe the cos 2(j) data, however, it fails 
to describe the function /i in a satisfactory manner. 

From the point of view of transverse momentum dependent distribution functions such a large cos 20 azimuthal 
dependence can arise at leading order only from a product of two T-odd functions, in particular, only from the 
distribution function h^ . 

We would like to mention the experimental observation that the cos 20 dependence as observed by the NAIO 
collaboration does not seem to show a strong dependence on A, i.e. there was no significant difference between 
the deuterium and tungsten targets. Hence, it is unlikely that it is in fact dominated by nuclear effects instead 
of effects associated purely to hadrons. Therefore, the unpolarized cross section as can be measured at RHIC 
is also likely to show a large cos 20 dependence, although replacing the pion by a proton will probably have a 
suppressing effect. 

Hence, we conclude that although there exist, apart from the h^ mechanism, several explanations of single 
spin asymmetries and also of the unpolarized cos 20 dependence in the Drell-Yan cross section, none of the 
approaches relate the two types of asymmetry and most of the effects are expected or found to be (too) small. 
Moreover, the effects should not only be large, they should also exhibit the right Qt behavior. 

III. TRANSVERSE MOMENTUM DEPENDENT DISTRIBUTION FUNCTIONS 

In this section we will discuss the transverse momentum dependent distribution functions that are needed to 
find the expressions for the leading order unpolarized and polarized Drell-Yan process cross sections differential 
in the transverse momentum of the lepton pair. 

We consider again Fig. ^. The momenta of the quarks, which annihilate into the photon with momentum 
q, are predominantly along the direction of the parent hadrons. One hadron momentum (Pi) is chosen to be 
along the lightlike direction given by the vector ri+ (apart from mass corrections). The second hadron with 
momentum P2 is predominantly in the n_ direction which satisfies n_)_ • n_ = 1, such that Pi ■ P2 = 0(q^). We 
make the following Sudakov decompositions: 



xiM^ 
5iV2'"+ ' QV2 



P|' = ^,,M+:iii:^„A^^ (6) 



"2 = 7=^ nl H -pz nz , (7) 

q^=^n^l + -^n^_+ ql^, (8) 



for Q^ = —q^ = q^ <^ Q^. We will often refer to the ± components of a momentum p, which are defined 
as p^ — p ■ nip. Furthermore, we decompose the parton momenta p, k and the spin vectors Si, S2 of the two 
hadrons as 

k^^n^ + '^'^^'+f^K ^ + kr, (10) 

X2V2 xQV2 

o - -^iQ xiXiMi 

oi = ^"-+ i^n^ + Sit, (H) 

X1M1V2 QV2 

„ _ A2Q X2X2M2 

The four-momentum conservation delta-function at the photon vertex is written as (neglecting 1/Q^ contri- 
butions) 



S\q -k~p) = S{q+ - p+) S{q- - /c") S^Pj. + fey ^ 9t): (13) 

fixing xP^ = p^ = q^ = xiP^, i.e., x = xi and similarly x — X2, and allows up to 1/Q^ corrections for 
integration over p^ and fc+. However, the transverse momentum integrations cannot be separated, unless one 
integrates over the transverse momentum of the photon or equivalently of the lepton pair. 

The parametrization of ^{p) should be consistent with requirements imposed on $ following from hermiticity, 
parity and time reversal invariance. The latter is normally taken to impose the following constraint on the 
correlation function p 14 



$*(P,5;p)=75C$(P,5;p)Ct75 (14) 

where p = (p", —p), etc. For the validity of Eq. (O) it is essential that the incoming hadron is a plane wave 
state. We will not apply this constraint and in the last section we will discuss this issue in detail. 

In the calculation in leading order we encounter the correlation function integrated over p~ , which is 
parametrized in terms of the transverse momentum dependent distribution functions as p7| ] 

dp' <i>(Fi, Si;p) = -^ { MxuPt) If + JM^uPt) e^.pal'^ 'Z^'^ 



-j^ - hM^uPr) j^^ h.^ixuPr) j^, + h, {xuPt) -j^ 



We used the shorthand notation 



ffis(a;i, Pt) =>^i9il{xi,Pt)-\ \f ^^ 9it{xuPt) (16) 



and similarly for hf-^. The parametrization contains two T-odd functions, that would vanish if the constraint 
Eq. (|lj) would be applied, i.e. the Sivers effect function f^ and the analog of the Collins effect, hj^. 
The parametrization of $ is 



Hx2,kT)^ I dk+ <^{P2,Sr,k) = ^ \f,{x2, fcT) ^ + /fT(^2, k^) e^.p.T^ ^^^^f ^ 

•' k~=X2P2 ,k^ ^^2 t ■^ ^ 

+9is{x2,krr) -j^ - hiT{x2,krr) ^^_^ /ii^(a;2, fcy) ^ ^^2 + ^1 i^2, k^) -^^jp — } ■ (17) 

The Sivers effect function f^ has the interpretation of the distribution of an unpolarized quark with nonzero 
transverse momentum inside a transversely polarized nucleon, while the function hj^ is interpreted as the 
distribution of a transversely polarized quark with nonzero transverse momentum inside an unpolarized hadron. 
In both cases the polarization is orthogonal to the transverse momentum of the quark. 

In terms of these functions we can schematically say that in order to fit the pp^ -^ n X data, Anselmino 
et al. |q| consider the following options for the product of three functions that are parametrizing the three 
soft parts: f^{xi,prp)fi{x2)Di{z) and hi{xi)fi{x2)H^{z,krp), where /i(x2) (or Di{z)) can be the gluon 
distribution (or fragmentation) function g{x2) (or G{z)) instead also. However, there is one remaining option 
(for pion production), that we are advocating as a source of single spin asymmetries: hj^{xi,pj,)hi{x2)Di{z). 
Because of the appearance of two chiral-odd quantities this contribution might be expected to be smaller 
than f^{xi,prp)fi{x2)Di{z). But even though /i > hi, one cannot exclude that hj^{xi,prp) is larger than 

Note that the magnitude of the Collins effect fragmentation function H^{z,krp) need not be related to 
the magnitude of hj^{z,kj,). In contrast to the distribution function, the fragmentation function will receive 
contributions due to the final state interactions which are present between the produced hadron and the other 
particles produced in the fragmenting of a quark. This is also the reason a similar constraint like Eq. ( [l4[) does 
not apply to fragmentation correlation functions. 



IV. UNPOLARIZED AND SINGLE SPIN DEPENDENT CROSS SECTIONS 

The DrcU-Yan cross section is obtained by contracting the lepton tensor with the hadron tensor 

U. +('^''')- (18) 



W^"" = ]- I dp-dk+cfprptfkrr S^{pt + krp- qj,) Tr ($(p) V^ $(fc) V^) 



The vertices Vl^ can be either the photon vertex V^ = e^^ or the Z-boson vertex V^ = gvl^ + 9Albl^ ■ The 
vector and axial-vector couphngs to the Z boson are given by: 

V, (19) 

(20) 



'3' 



9a — ^3 



where Q^ denotes the charge and T| the weak isospin of particle j (i.e., T| = +1/2 for j = u and T| — —1/2 
for j = e^,d,s). We find for the leading order unpolarized Drell-Yan cross section, taking into account both 
photon and Z-boson contributions. 



da^°Hhih2 -> eix) _ oP- 

d^dx\dx2d^ qrj^ 3Q^ 



Y^ A'i(2/).F[/iA] 



+ [X3(y)cos(2,^) + X4(y)sin(20)] T 

and for the case where hadron one is polarized: 

da^^'>{hih2^£eX) _ a2 r 

dfldxidx2d^qj, 3Q^ ^ | 
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-\Sit\ [-f'^3(y)sin(30-(/)sJ -i^4(y)cos(3(/«-(?!)Si)] 
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Ah 



Ah-kj.{h-pTf -2h-prj, pj,-kj. - hkj. p| j 2M'^M 



(22) 



where the ellipsis stand for the T-even-T-even structures (which for the contribution of the virtual photon are 
absent, cf. Ref. jlj]). Let us list the various definitions appearing in these expressions. We have defined the 
following combinations of the couplings and Z-boson propagators: 



Ki{y) = A{y) [el + 2g\,eag'^Xi + 4c?X2] 



c{v) 



[2g\eag'\xi + 4^3X2] 



K2{y) = A{y) [2g\,eag'Xxi + ^cgxa] 



C{v) 



\2g\eagvX\ + c\c\x'^ , 



Kz[v) = B{v) {el + 2gi.ea5yXi + ^0^X2] 
K^iy) = B{y) {2g'yeag\x^ , 
which contain the combinations of the couplings 

4 = ^SvS'a- 



d^{gi'-g\'], 



j — i or a 



(23) 

(24) 

(25) 
(26) 

(27) 
(28) 
(29) 



The Z-boson propagator factors are given by 

1 



i2/,n2 
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X3 



(Q'-M^,) 



sm'{2ew) (Q2-Af|)2 + r|Af|' 
1 Q" 
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(30) 
(31) 
(32) 



The above is expressed in the so-called CoUins-Soper frame |28|, for which we chose the following sets of 
normalized vectors (for details see e.g. [pj[): 



i = q/Q, 

h = qr/Qr = {q - Xl Pi ~ X2 P2)/Qt, 
where Pi = Pi — q/{2xi), such that: 



1 
1 

71 



f + zf" 
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(33) 
(34) 

(35) 

(36) 

(37) 



The azimuthal angles lie inside the plane orthogonal to t and z. In particular, dfl — 2dyd(j)'', where (j)^ gives 
the orientation of l'^ = [g^'^ — t^^t^ + z'^z'^) l,y, the perpendicular part of the lepton momentum I; (p, (ps^ are the 
angles between h,S^rp and l±, respectively. In the cross sections we also encounter the following functions of 
y = l^ /q^, which in the lepton center of mass frame equals y = {1 + cos6')/2, where 9 is the angle of z with 
respect to the momentum of the outgoing lepton I (cf. Fig. g): 



A(2/)=Q-y + 2;y=i(l + cos2. 
i?(y) = y(l-Z/) = isin^e, 



C{y) = {l-2y)'^'-cos0. 
Furthermore, we use the convolution notation (Ralston and Soper [pi use /[•..]) 



J^lff] = J d'pT d'kj. S'ipr + kr - qT)nxi,pj^)nx2,k'T), 



(38) 
(39) 
(40) 

(41) 



where a is the flavor index. 




lepton plane (cm) 
FIG. 3. Kinematics of the Drell-Yan process in the lepton center of mass frame. 
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Since we are mainly concerned with the single polarized Drell-Yan process, we have given the double polarized 
cross section in the appendix for completeness and future reference. 



V. QUARK SPIN CORRELATIONS 



In order to explain the angular dependence of the unpolarized cross section as measured for the process 
7r~ TV ^ /i+ /i~ X , where N is either deuterium or tungsten, using a 7r~-beam with energy of 140, 194 and 286 
GeV H] (—1/2 < cos6' < 1/2), Brandenburg et al. [g^ proposed factorization breaking correlations between the 
transverse momenta of the incoming quarks and between their transverse spins. This correlation between the 
transverse momenta is taken to be 



P{Px, krp) (fpj, (fkj. 



aT{oLT + 2/3^) 



exp [-aripT + ^t) ^ Pt{Pt ~ ^tT] d^Pr cPkj., 



(42) 



which reduces to separate Gaussian transverse momentum dependences, since j3t is found to be practically zero 
(and ay = 1 GeV~^ at these energies). 

In case the boson V that produces the lepton pair, is a virtual photon {V = 7*) Brandenburg et al. fit the cross 
section Eq. (pi) using the data for the 7r~-beam with energy of 194 GeV and lepton pair mass m-y* = 8 GeV/c^. 
They find that 1 — A — 2z/ « —4k, where they take the following model for k, which is a measure of the correlation 
between the transverse spins of the incoming quarks, 



Ko 



Qt 



Q: 



(43) 



and the fitted values are kq = 0.17 and rriT = 1-5 GeV. 

We will do a similar analysis based on the assumed presence of T-odd distribution functions with intrinsic 
transverse momentum dependence. For simplicity we take /i = 0, A = 1 (in accordance with the expectation 
from next-to-leading order perturbative QCD and the data in the Collins-Soper frame) and define i^ = 2k. For 
y = 7* we then find the following expression for n (cf. Eqs. (h) and (Ell) 



E^a 



T 



A ft ' P'T' ft ' fZrp P'T' ' "'7 



h^hf 
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E-^^^tA/i] 



(44) 



A model for the shape of the function hi is needed. Collins' parametrization JTJ for the fragmentation function 



Hi is (note that Collins uses the function AD 



H/a 



e^SiTikTjH{ 









Ml 



Im[A*(fc2)B(fc2 



(1-z) 



(45) 



where Mh is the mass of the produced hadron and in his model Mc is the quark mass that appears in a dressed 



fermion propagator i{A{k^)^ -I- B{k^)Mc)/{k^ 



Ml) 



the functions A and B are unity at k = M^ 



We assume a similar form for hi in terms of /i (we assume no flavor dependence of Mc)'. 



/ijL'^ 



{x,Pt) 



McMh 



fHx^p'r) 



^"W 



Ml" 



(46) 



using the constant djj (which in principle is a function of x) and Mc as the fitting parameters now (and similarly 
for the antiquark distribution functions). We also assume the above given Gaussian transverse momentum 
dependence for /i(a:,p|n). After multiplying Eq. (H) by a trivial factor Q\IQ\^ using the fcr integration to 
eliminate the delta function and shifting the integration variable pt -^ p'x — Pt ~ ^Qti ^^^ arrives at 



Ki ax 



IT 



MlQ 



j d^PT 



{p't + Ht) 



Ml 



(p't - Ht)' 



Ml 



-2aTPT 



(47) 



where ki = c/f^ 0^^/2/2 and for the moment we considered the one flavor case. We approximate this by taking 
p't ~ (where the exponential factor is largest) in the term between square brackets (this is valid for large 
enough values of Qt, but the resulting expression also has the right Q^ behavior as Qt — > 0), this results in 
(reinstalling the flavor summation) 



Q|M^ 



Ea,aeg^?/r(^o/?(^2) 

Let us for simplicity also assume kJ to be independent of the flavor and fit 



8ki 



QlAPc 



[Ql 



-\MlY 



(48) 



(49) 



to the data at 194 GeV of Ref. [||. This does not give as good a fit (Fig. Q) as a factor of Q\ in the numerator 
would give (for this particular set of data), but it can obviously reproduce the tendency. Moreover, it has the 
desired property that k vanishes in the limit of Qt ^ oo, as opposed to Eq. (43). We find for the dressed quark 
mass Mc a rather large value of 2.3 ± 0.5 GeV compared to the chiral symmetry breaking scale, but one should 
not take the model too seriously and we have made several approximations. 




Qt [C5eV] 

FIG. 4. Data from |3[ at 194 GeV and fit (using Eq. (M)) to t' = 2k as a function of the transverse momentum Qt 
of the lepton pair. The fitted parameters are Mc = 2.3 ± 0.5 GeV and 16 ki = 7 ± 2. 

We have chosen the data at 194 GeV of Ref. [g|, because it has the smallest errors (the error in Qt is chosen 
to be the bin size). The fits to the three other available sets of data, namely at 140 and 286 GeV of Ref. H 
and at 252 GeV of Ref. [Q, yield lower values of Mq and n\ (on average a factor 2 smaller), hence, have a lower 
maximum (at a smaller value of Qt) f^nd are less broad. We take the above result as providing a rough upper 
bound. 

Taking for simplicity c%. — c"^ — c°j^ , we arrive at a (crude) model for the function h\{x\^'p^\. 



ht 



aT a McMh 



{x,Pt) — Cfj 2 



M^ 



-o-tPt 



Mx), 



(50) 



with Ale — 2.3 GeV, c'^ — I and ar = 1 GeV^^, which can be used to get rough estimates for other asymmetries. 
The factor aTJT^ comes from the consistency requirement between the definitions of f\(x) and J\{x^ fc^^) with a 
Gaussian k\ dependence. In the next section we will discuss the relevant asymmetries for RHIC. 



VI. IMPLICATIONS FOR RHIC 



From Eq. (|2^) we see that in case V = ^* and that when we neglect the "higher harmonic" term containing 
the 3(^ dependence, there are two single transverse spin azimuthal dependences, namely sin((/) — 4>S\') arising 
with the Sivers function /j^ and sin((/) + c^Sx ) arising with h\; . 
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To estimate the size of the sm{(j) — (ps^ ) term, one can use for instance one of the usual parametrizations of 
/i and the parametrization for f^ as found by Q. To estimate the size of the sin(0 + (f>Si) term, one can use 
one of the models for hi [ [l2|p9| ] or take the upper bound for hi that arises from Soffer's inequality {gi is also 
well-known) and one can use a fit for h^ from the unpolarized azimuthal cos 2(j) dependence of the cross section 
in pp — > ££X, in a similar way as was done in the previous section. 

Let us examine the sin(0 + (j>Si) dependence of the cross section with the above given model for hj^. The 
relevant expression for the cross section in the polarized case is given by (cf. Eq. (^ with ^ = and A = 1) 

— —-^ — ex (l + cos^9 + K sin^9cos2(f)- p |S,y| sin^ 6* sin(0 + 05j + ...), (51) 

a ail dms. 



where the ellipsis stand for the other angular dependences. The analyzing power p is found to be (cf. Eq. (|22 



a. a 



ft ' fin~' 



M2 

Using the above model Eq. (M) for hj^ and performing similar approximations as before, we arrive at: 



Qi+AM'C Earaeim^l)ni^2) 2V«max Ea . ^^ /f (^l) /? (^2) 



(53) 



where Kmax is the maximum value of k, which is at Qt = 2Mc- A determination of p, k, hi should be mutually 
consistent according to the above equation, if the underlying mechanism is indeed the one that is assumed 
here. The maximum value p is also at Qt = 2Mc, which in the case of one flavor corresponds to Pmax = 
ch hi{xi) / {2fi(xi)) < ch/2 sa 1/2. If hi is for instance an order of magnitude smaller than /i, this would give 
an analyzing power for this single transverse spin azimuthal asymmetry at the percent level. 

The above scheme entails many extrapolations and assumptions and prevents us from stating accurate esti- 
mates for the asymmetries. One problem comes from the fact that the fit in the previous section resulted from 
data of the process 7r~ A^ — > /x+ /^~ X, where N is either deuterium or tungsten, so extrapolation to pp ^ £iX 
is unclear. One might expect that the cos2(/) dependence of pp ^ i£X a,s will be measured at RHIC is smaller 
than for the process tt" N — ^ /i"*" p~ X, since in the former there are no valence antiquarks present. In this 
sense, the cleanest extraction of hj^ would be from pp ^ i£X. 

Another problem concerns the energy scale. The extrapolations should involve evolving the functions to the 
relevant energies, however, the evolution equations for f^ and h^ are not yet known. 

However, the basic idea is clear. One fits the unpolarized azimuthal cos 2(f) dependence of the cross section in 
a similar way as was done above, (for instance) by using a Collins' type of Ansatz to arrive at a model for hj^, 
which then can be used to measure or cross-check the function hi by measuring the sin((/) + (psi ) dependence. 

VII. DISCUSSION 

A cos2(j) term in the hadron tensor is itself a T-even quantity, but in our approach it is factorized into a 
product of two T-odd functions. From the definition of the correlation function $(p) one can show that time 
reversal symmetry requires the T-odd functions to be zero M. This assumes that the incoming hadrons can 
be describe as plane waves states. To circumvent this conclusion one could think of initial state interactions 
between the two incoming hadrons [pj or one could think of effects due to the finite size of a hadron pO] . 

Initial state interactions between the two incoming hadrons would be a factorization breaking effect (not to 
be confused with the breakdown of factorization at higher twist [ pl| ) and this implies nonuniversality of the 
functions involved. The factorization breaking correlations proposed by Brandenburg et al. 123] assuming some 
nonperturbative gluonic background |p2| , might be universal in some restricted sense. For instance, one could 
retain universality among a subset of possible processes, namely the ones with exactly the same initial states. 
This would mean that functions obtained from the Drell-Yan process can be used to predict asymmetries in 
the process pp — > ttX successfully. Another type of universality would be that the factorization breaking 
correlations are the same for different asymmetries in the same process, e.g. the same for v and p in the case 
discussed above. These issues can be tested experimentally. We have proposed a concrete way to test some of 
these issues. 

At finite scales Qt and Q one expects the finite size of a hadron to play a role. However, such nonperturbative 
effects should not conflict with the factorization formula for Drell-Yan at finite Qt and Q {Qt ^ Q) Q- The 
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finite size of hadrons most likely results in higher twist contributions, but maybe it will just prevent the naive 
application of Eq. ( |l4|) as a constraint imposed by time reversal symmetry, which would not conflict with the 
factorization formula. These issues need to be investigated further theoretically. 

Let us just mention that finite size effects have been proposed as origins for the Sivers effect in Ref. |30|. 
Spin-isospin interactions have also been proposed pj] to obtain a nonzero Sivers function. Liang et al. [ p5[ have 
proposed a model relating the spin of a hadron to the orbital motion of quarks inside that hadron. This could 
be viewed as a model for the function f^ and a similar model might be constructed for the function h^ . 

It is worth emphasizing that the functions f^ and h^ appear in quite different asymmetries in general, even 
though they can both account for the single spin asymmetries in pp^ -^ tt X. For instance, f^ cannot account 
for the cos 20 asymmetry discussed above and also, it yields a different angular dependence for the single spin 
asymmetry in the Drell-Yan cross section as was pointed out in the previous section. Also, in contrast to fef , 
the Sivers effect, which is chiral-even, might produce single spin asymmetries in (almost) inclusive DIS ||6|, p4| , p6[ , 
unless it originates from initial state interactions between hadrons. 

It is good to point out that the Berger-Brodsky higher twist mechanism is not ruled out as a possible 
explanation for the observable i', although in the higher twist model of Ref. |23tl using a pion distribution 
amplitude it seems to fall short in explaining the large values found for v. Of course, it might contribute in 
addition to the hj^ mechanism. However, an observed correlation between v and p will be indicative of the 
latter. 



VIII. CONCLUSIONS 

We have discussed in detail the consequences of T-odd distribution functions with intrinsic transverse mo- 
mentum dependence for the Drell-Yan process. In particular, we focused on a chiral-odd T-odd distribution 
function, denoted by hj^, which despite its conceptual problems, can in principle account for single spin asym- 
metries in pp^ — > ttX and the Drell-Yan process, and at the same time for the large cos 20 asymmetry in the 
unpolarized Drell-Yan cross section as found in Refs. P,0, which still lacks understanding. We have used the 
latter data to arrive at a crude model for this function and have shown explicitly how it relates unpolarized 
and polarized observables that could be studied at RHIC using polarized proton-proton collisions. It would also 
provide an alternative method of gaining information on the transversity distribution function hi. 

The distribution function hj^ would signal a correlation between the transverse spin and the transverse 
momentum of quarks inside an unpolarized hadron. It is formally the distribution function analog of the Collins 
effect, which concerns fragmentation, but most likely arises from quite a different physical origin. Further 
theoretical and experimental study of these issues is required. 

We have also listed the complete set of azimuthal asymmetries in the unpolarized and polarized Drell-Yan 
process at leading order involving T-odd distribution functions with intrinsic transverse momentum dependence. 
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APPENDIX 

The leading order double polarized Drell-Yan cross section, taking into account both photon and Z-boson 
contributions, is found to be 



da^^\hih2 -^ UX) 



a 



2 



a, a 



.^^ \Sij.\\S,j.\ cos(20~05i"0sj^ 
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(54) 



where the ellipsis stand for the (remaining) T-even-T-even structures (which for the contribution of the virtual 
photon can be found in Ref. |l4|). 
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